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Abstract. We analyze the behavior of the eigenvalues and eigenfunc- 
tions of the Laplace operator with homogeneous Neumann boundary con- 
ditions when the domain is perturbed. We show that if f2o C f2e are 
bounded domains (although not necessarily uniformly bounded) and we 
know that the eigenvalues and eigenfunctions with Neumann boundary 
condition in Q e converge to the ones in Qq, then necessarily we have that 

|f2 E \ f!o| — > while it is not necessarily true that dist(f2 e , Qq) 0. As 
a matter of fact we will construct an example of a perturbation where 

the spectra behave continuously but dist(O e ,f2o) +co. 



1. Introduction 

This paper is concerned with the behavior of the eigenvalues and eigenfunc- 
tions of the Laplace operator in bounded domains when the domain undergoes 
a perturbation. It is well known that if the boundary condition that we are 
imposing is of Dirichlet type, the kind of perturbations that we may allow in 
order to obtain the continuity of the spectra is much broader than in the case of 
Neumann boundary condition. This is explicitly stated in the pioneer work of 
Courant and Hilbert [5] and it has been subsequently clarified in many works, 
see [U [2j [6] and reference therein among others. See also [8] for a general text 
on different properties of eigenvalues and [S] for a study on the behavior of 
eigenvalues and in general partial differential equations when the domain is 
perturbed. 

In particular, with Dirichlet boundary condition we may consider the case 
where the fixed domain is a bounded "smooth" domain f2o C M , N > 2, and 
the perturbed domain is f2 e in such a way that Qq C f2 e , that is we consider 
exterior perturbation of the domain. We may have perturbations of this type 
where |f2 e \ Qo| > V for some fixed r/ > and still we have the convergence 
of the eigenvalues and eigenfunctions. Moreover, we may even have the case 
|O e \ f2o| — > +oo and still we have the convergence of the eigenvalues and 
eigenfunctions. 
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To obtain and example of this situation is not too difficult. If we consider 
for instance fi C M 2 , given by fio = (0, 1) x (—1, 0) and 

fi e (a) = {(x, y) : < x < 1, -1 < y < a(l + sin(x/e))} D fi 

where a > is fixed, we can easily see that the eigenvalues and eigenf unctions 
of the Laplace operator with Dirichlet boundary condition in fi e converge to the 

ones in fio. Moreover |fi e | = |fio| + J Q a(l + sin(x/e))dx ~ |fio| + a for e small 
enough. Moreover, it is not difficult to modify the example above choosing the 
constant a dependent with respect to e in such a way that a(e) — > +oo and such 
that we still get that the eigenvalues and eigenfunctions in fi c (a(e)) converge 
to the ones in fi an d |fi c (a(e)) \ fio| — ► +oo. This example shows that the 
class of perturbations that we may allow to get the "spectral convergence" of 
the Dirichlet Laplacian is very broad and that knowing that the eigenvalues 
and eigenfunctions of the Dirichlet Laplacian converge does not have many 
"geometrical" restrictions for the domains. 

The case of Neumann boundary condition is much more subtle. As a matter 
of fact, for the situation depicted above it is not true that the spectra converge. 
So we ask ourselves the following questions: if we have a domain fio and 
consider a perturbation of it given by fio C fi e , where we assume that all the 
domains are smooth and bounded although not necessarily uniformly bounded 
on the parameter e, then if we have the convergence of the eigenvalues and 
eigenfunctions, 

(Ql) should it be true that |fi e \ fi | ^ 0? 

(Q2) should it be true that dist(fi e ,fio) = sup xef2 dist(x,fio) — 0? 

We will see that the answer to the first question is Yes and, surprisingly, 
the answer to the second one is No. 

Observe that, as the example above shows, the answer to both questions for 
the case of Dirichlet boundary condition is No. 

In Section [2] we recall a result from (TJ [3] which provides a necessary and 
sufficient condition for the convergence of eigenvalues and eigenfunctions when 
the domain is perturbed. In Section[3]wc provide an answer to question (Ql) 
and in Section 3] we provide an answer to question (Q2). 

2. Characterization of spectral convergence of Neumann 

Laplacian 

In this section we give a necessary and sufficient condition for the con- 
vergence of the eigenvalues and eigenfunctions of the Laplace operator with 
Neumann boundary conditions. We refer to [T] and [3] for a general result in 
this direction, even in a more general context than the one in this note. In 
our particular case, we will consider the following situation: let fio be a fixed 
bounded smooth (Lipschitz is enough) open set in M. N with N > 2 and let fi e 
be a family of domains such that for each fixed < e < eo, fi e is bounded and 
smooth with fio C fi e . 

Let us define now what we mean by the spectral convergence. For < 
e < eo, we denote by {A^ 1 }^ 1 the sequence of eigenvalues of the Neumann 
Laplacian in fi e , always ordered and counting its multiplicity, and we denote 
by a corresponding set of orthonormal eigenfunctions in fi e . Also, 
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since we are considering domains which vary with the parameter e and we will 
need to compare functions defined in fJo and in O e , we introduce the following 
space Hi = H 1 ^) © H 1 ^, \ O ), that is x G H\ if X\Q G ^(^o) and 
X|(n e \n ) e ^(^e \ Oo), with the norm 

11x11*} = llx||iri(no) + llx||^i ( n e \fio) ' 
We have that iJ 1 (S7 e ) °-> if* and in a natural way we have that if X G 
i/ 1 (fio) via the extension by zero outside fio we have x G i?*. Hence, with 
certain abuse of notation we may say that if Xe G if,}, < e < eoj then 

Xe ^ Xo in if ||xe - Xo\\m(n„) + ||Xe||ffi(n«\fio) ^ °- 



Definition 2.1. FKe will say that the family of domains Q e converges spectrally 
to fio as e — > i/ i/ie eigenvalues and eigenprojectors of the Neumann Laplacian 
behave continuously at e = 0. That is, for any fixed n G N we have that 
A* — > A° as e — - > 0, and /or eac/i n G N swc/i that A„ < A° +1 £/ie spectral 
projections : L 2 (M JV ) -» tf 1 ^), P^(V>) = E2=i V0/i'(n.)#, satisfy 

S up{||P^) - P n °(V)ll^, ^ G L 2 ^), = 1} ^2 0. 



The convergence of the spectral projections is equivalent to the following: 
for each sequence ek —> there exists a subsequence, that we denote again 
by e^. and a complete system of orthonormal eigenfunctions of the limiting 
problem such that ||<^ - \\ H i -> as k -> oo. 

In order to write down the characterization, we need to consider the follow- 
ing quantity 

' |V0| 2 



min . (2.1) 



6=0 111 Q 



Observe that r £ is the first eigenvalue of the following problem with a combi- 
nation of Dirichlet and Neumann boundary conditions: 

-Aw = tu , Q e \ Clo , 
u = , <9f2o , 

— = o , an £ \ 5fi • 

an 

We can prove the following, 

Proposition 2.2. ^4 necessary and sufficient condition for the spectral con- 
vergence of f2 e to f2o * s 

r e ^+oo. (2.2) 

We refer to pQ and [3] for a proof of this result. 

Remark 2.3. The fact that Qq C f2 c can be relaxed. It is enough asking 
that for each compact set K C Oo there exists e(K) such that K C fi e /or 
< e < e(#), see [3]. 
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3. Measure convergence of the domains 

In this section we provide an answer to the first question. Observe that in 

Proposition 1 2 . 21 we do not require that \Q e \ Q \ ^— S 0. However, we have the 
following 

Corollary 3.1. In the situation above if fi e converges spectrally to Qq, then 
necessarily |O e \ fiol 0. 

Proof. This result is proved in [5] but for the sake of completeness and since 
it is a simple proof, we include it in here. 

If this were not true then we will have a positive 77 > and a sequence Ek — > 
such that |r2 £fc \ f2o| > ?7- Let p = p(r)) be a small number such that \{x G 
R N \O , dist(x,fi ) < p}\ < V/ 2 - Tnis implies that \{x G Q Ck , dist(x,O ) > 
p}\ > ? //2- Let us construct a smooth function 7 with 7 = in f2o, and 
7(x) = 1 for x G R N \ flo with dist(a;,Oo) > P- Then obviously 7 G H 1 ^^) 
with ||V7||i2 (nefc ) < C and ||7||i3(n eh ) > (v/ 2 )^- This implies that r £fc is 

bounded. Hence it is not true that r e ^5 +00 and therefore, from Proposition 
12.21 we do not obtain the spectral convergence. ■ 

In particular, this result implies that the answer to question (Ql) is af- 
firmative. That is, if we have the convergence of Neumann eigenvalues and 

eigenfunctions, necessarily we have that |Sl e \ Oo| > 0. 

4. Distance convergence of the domains 

In this section we will provide an answer to question (Q2) and, as a matter 
of fact, we will see that the answer is No. We will prove this by constructing an 
example of a fixed domain f^o and a sequence of domains f2 e with Qq C fl e with 
the property that dist(f2 e , flo) does not converges to 0, but the eigenvalues and 
eigenfunctions of the Laplace operator with Neumann boundary conditions in 
tt e converge to the ones in fig, see Definition 12. II 

As a matter of fact in [31 Section 5.2] a very particular example of a dumbbell 
domain (two disconnected domains joined by a thin channel) is provided so 
that the eigenvalues from the dumbbell converge to the eigenvalues of the 
two disconnected domains and no spectral contribution from the channel is 
observed. In this note we will obtain a family of channels for which the same 
phenomena occurs, see Corollary 14.41 and will provide a proof, different from 
the one given in [3J . 

Let us consider a fixed domain Slo C t 1 *' which satisfies that ilo c {x G 
R N , xi < 0} and such that 

O n {x = (x u x') G K x IR^ -1 , -1 < x x < 1, \x'\ < p} 

= {x= (xi,x') G K x l" -1 , -1 < xi < 0, \x'\ < p} 

for some fixed p > 0. 

We will construct £l e as f2 c =int(Oo U -R e ), where R e is given as follows 

R e = {(x u x') elx : < x x < L, \x'\ < g e (xx)} (4.1) 
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where the function g e will be chosen so that g e > 0, g e £ C 1 ([0, L]) and 
g e — > uniformly on [0, L], sec Figure [T] For the sake of notation we denote 
by T e = dR e n {x x = 0} and T e L = dR e n {x l = L}. 




Figure 1. The exterior perturbation R e . The thick line refers 
to the supplementary Dirichlet condition in the problem (|4.2[) , 
while Neumann boundary conditions are imposed elsewhere. 



We refer to [T^] for a general reference on the behavior of solutions of partial 
differential equations on thin domains. See also the recent survey [7] for a 
study on the spectrum of the Laplacian on thin tubes in various settings, and 
for many related references. 

Observe that if L is fixed then dist(Sl e , Slo) = L for each < e < eq. 
Moreover we will show that for certain choices of g e we obtain the spectral 
convergence of the Laplace operator. To prove this results, we use Proposition 
2.21 and show that r e — > +oo. Notice that r e , defined in (|2.ip is the first 
eigenvalue of 

-Ait = tu, R t , 

u = o, r§, £ 4>2 j 

| = o, d R A n. 

Since we have Neumann boundary conditions on the lateral boundary of i? e , 
there clearly exist profiles of g e for which r e remains uniformly bounded as 
e — > 0. In fact, a simple trial-function argument shows that r c < ir 2 /(2L) 2 
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whenever g e (s) > g e (Q) for every s G [0, L\. The idea to get r e — > +00 consists 
in choosing a rapidly decreasing function s 1— > g e (s), which enables one to 
get a large contribution to r e coming from the longitudinal energy due to the 
approaching Dirichlct and Neumann boundary conditions in the limit e — > 0. 
Let us notice that a similar trick to employ the repulsive contribution of such 
a combination of the boundary conditions have been used recently in [TO] to 
establish a Hardy- type inequality in a waveguide; see also for eigenvalue 
asymptotics in narrow curved strips with combined Dirichlct and Neumann 
boundary conditions. In our case, we are able to show 

Proposition 4.1. With the notations above, for any function 7 £ C 2 ([0, L]) 
satisfying 

< a < 7 < ai < 1, j(L) < 0, and 7 > a 2 > (4.3) 

for some positive numbers ao, ct\ and a 2 , if we define g t — j 1 ^ we have that 
r e — > 0. 

In particular, applying Proposition \2.2\ we obtain the convergence of the 
eigenvalues and eigenf unctions of the Neumann Laplacian in fl e to the ones in 
O . 

Remark 4.2. Observe that a function 7 satisfying (|4.3[) necessarily satisfies 
that j(s) < for < s < L. Hence, the function 7 is decreasing. 

Proof: Since r e is given by minimization of the Rayleigh quotient, 

T e = inf 

4>€H 1 (R 

0=0 in r, 







we analyze the integral J R |V(/>| 2 for a smooth real- valued function (f> with 
<f> = in a neighborhood of T^. We have 

/ |V0| 2 =/ f {\<f> Xl \ 2 + \V x ,<t>\ 2 )dx'd Xl 

J R, Jo J\x'\<g € (x 1 _) 

Considering the change of variables x\ = y\, x' = ge{yi)y' which transforms 
(xi,x') s R € into (yi,y') G Q where Q is the cylinder Q = {{yi,y') : < 
yi < L, \y'\ < 1} and performing this change of variables in the integral above, 
elementary calculations show that 

2 



f f ( ' N \ 1 N 



9 i r i dy 



where </j(y) = <f>(yi,g e (yi)y')- 

N — 1 

Writing the above expression in terms of the new function tp(y) = ge{yi)~ ~ f(y) 
so that 

9l N - l)/2 <Pvi=4>vi> i = 2 t ...,N, 
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we get, 



IW| S 



> 



Q 



N - 1 ' ' N \ 2 1 N 

Z 5c 5e J= 2 / ffe *=2 



2 5e 



V 



3e 



i=2 
N 



■ 2 N N 



& ,=2 
2 A.2 



i=2 



-V - (iv- i)— Ma 

9e 



97 



N 



N 



+ ( N - 1 )^Y,y^ + ij2< d y 

9e i=2 yc i=2 

where we have used that (ip Vl — 2i=2 Vi^Vi^) 2 — 0- Via integration by parts 
in the second and third term above, we get, 

r-L 



y'\<l JO 



-(N-l)^(ip 2 ) yi dyidy' 



\v'\<i 



yi=L 



yi =o JO 



y'l<i 



(N-l)M^-^(L,y')dy' + 



.5e 5e 



and 



/• -2 w r L a 2 N r \ 

/ (iV - 1)% V y^ v Mv = (N - 1)% J2 / Vioi^dl/dy 
JQ 9e l=2 Jo 9e J^J\y'\<i z 



N-l g 2 
2 5? 



y' l=i 



■0 — (iV — 1) / dj/i 



Iv'Ki 



Hence if we require that ge(L) < 0, we have, 



|V«^| 2 > 



Q 



N-l g t 

2 ffe 



2 5? 



iV - 1 V iV - 1 \ gl 

N 



ip 2 dy' I dyi + 



|y'l=i 



JQ 9e i=2 



(4.4) 
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The last two terms in this expression can be written as 



1 



and we have that 



l\y'\<l 



N - 1 



N - 1 



■ 9e 



y'l<i 



9l(Vi) 



^ 2 dyi 



li/'l=i 



y'\ 



1? 



|y'l<i 



with p = p(yi) being the first eigenvalue of the problem 

-Ay.]) = fn/> , \y'\<l 




i 



where n denotes the outward unit normal vector field to the (N—2) dimensional 
unit sphere Si = {y' G : \y'\ = 1}. 

We claim that if we denote by A (77) the first eigenvalue of 

-A v 4 = \il>, \y'\<i, 



■5- + w 

on 



0. 



1 



we have that 



1 5 1 I 

— r j^Yj as 77 — > 0, where Bi is the (N — 1) dimensional unit 

ball and Si its surface, which satisfy |5i| = (AT — l)|v3i|. As a matter of fact 
by standard continuity result we know that A (77) — * and its eigenfunction 
t/>, ( , which is radially symmetric, converges to the constant function 1/ ^J\B\\, 
which is the first eigenfunction of the Neumann eigenvalue problem. But 



which implies that 



Bi 



iv^l 



> 



Si 



\Si 



Si 



Moreover, using ijj = l/y/\Bi\ as a test function in the Raylcigh quotient 

This proves our claim. In 
= r]o(5) such that A (77) > 



. \Si\ 

\Bi 



for A(?7), we immediately obtain A(?7) < 77 

particular, given 6 > small, we can choose 770 
(N - 1 - S)n for < 77 < ?i . 

Therefore, if we choose the function g e such that g e (yi) 



uniformly in 



yi G [0, L], we have that p(yi) > ^ N 1 ^ Q ^ivi) f° r 6 small enough. 



Hence, 



|V<^| 2 > 



N — 1 



2 <? e 

(AT- 1)(AT- 1 -6) N-l 



ip dy 



N - 1 



AT- 1 



(N — 1 — 6) + 1 
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and observe that the number k = ^— ^ — (iV — 1 — 5) + 1 is strictly less than 
one for all values of N > 2 choosing a fixed and small <5 > 0. If we denote by 



inf (* 

0<yi<i \g e 



then 



|V0| 2 > 



N - 1 



AT- 1 



N-l 



Consequently, r e > 

Let us see that we can make a choice of the family of functions g e , satisfying 
the two previous conditions we have imposed, that is g t {L) < and g e (vi) — * 
uniformly in < y\ < L such that m e — > +oo as e — > 0. 



Let us choose a function 7 £ C 2 ([0, L\) satisfying 
Then, we have 



and let g e 



7 



1 1-1. 

-7« 7: 



9t 



1/1 s 1 9 1 

_(__l) 7 7-2^ 



- 7 e 7) 



and simple calculations show that 



fj< 



1,1 



-(--!)-« - 



67 ~ a e 



for e > small enough so that — 1) 
m e — > +00 and it proves the proposition. 



> . This shows that 



Remark 4.3. Now that we have been able to construct a thin domain R £ as 
in (|4.1[) such that T e ^5 +oo ; we can construct another thin domain R e such 
that its "length" goes to infinity, its width goes to zero and still f t ^5 +00, 
where f e is the first eigenvalue of (|4.2p in R t instead of R e . 

For this, let R e be a thin domain constructed as in Proposition \4-l\ and let 
p e be a sequence with p e — > +00 such that — > +00 and ct 1 ^ p e — > 0. Define 

R e = p e R e , that is 

R f _ = {{xi,x') : < x\ < p e L, \x'\ < p t g e (xi)} , 

then < p e g e (x\) < a\ ' e p e < -^> and f e = ?± +00. 

Observe that if we require also a Dirichlct boundary condition in T e L , we 
can relax the conditions on 7 in Proposition ^. H and in particular the condition 
i{L) < can be dropped. Hence, we can show, 

Corollary 4.4. With the notations above, for any function 7 £ C 2 ([0, L]) 
satisfying 

< cto < 7 < ai < 1, and 7 > c*2 > 
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for some positive numbers olq, otx and a 2 , if we define g ( = j 1 ^ we have that 
f e — ^> 0, where ft is the first eigenvalue of 

Re , 

Proof: This follows easily by a Neumann bracketing argument. More pre- 
cisely, from the hypotheses, 7 is a strictly increasing function. Hence, cither 
7 is strictly monotone in (0, L), or there exists a unique L* € (0, L) such that 
7 (£*)=0. 

In the first case, if 7 is decreasing (respectively increasing) we substitute the 
Dirichlct boundary condition at (respectively at Tq) by a Neumann one. 
Then the new eigenvalue problem gives rise to r e defined exactly in the same 
way as (|4.2[) (modulo possibly a mirroring of R e ) and we have f e > r e — » +00 
as e — > 0. 

In the second case, we cut the domain R e in two domains R® = R e n {0 < 
x x < L*}, R\ = R c (7 {L* < xi < L}. We know that f e > inf{r e °, r^}, where 
r e ° and r e are the corresponding eigenvalues in i?^ and iZ^ with a Neumann 
boundary condition imposed at the newly created boundary R e n {x\ = L*} 
on both domains. In both domains we can apply Proposition l4.1l as in the first 

case so that r° , t\ ^> +00, which implies f e — > 0. ■ 

Remark 4.5. XTus corollary recovers and generalizes the results from Sec- 
tion 5.2 in [3]- 
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